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Research  on  this  project  has  centered  on  the  Phi-transform  and  wavelets  and  their 
application  to  various  problems  in  image/signal  processing,  differential  equations,  and 
computer  aided  design. 


Construction  of  wavelets.  While  the  Phi-transform  and  wavelets  have  been  developed 
in  a  quite  general  setting,  several  applications,  most  notably  in  image  processing  and 
numerical  methods  for  differential  equations,  have  pointed  to  several  deficiencies  in  the 
traditional  wavelet  constructions.  This  has  led  us  to  consider  alternate  constructions  of 
wavlelets  and  related  questions. 

In  a  series  of  papers  [6,7,8],  C.  de  Boor,  R.  DeVore  and  A.  Ron  study  the  structure  of  shift 
invariant  spaces  which  play  a  prominent  role  in  wavelet  constructions,  multigrid  methods, 
applications  of  radial  basis  functions,  and  various  aspects  of  approximation  (splices  of 
univariate  and  multivariate  splines  often  have  this  property).  A  space  S  of  functions  on 
Rd  is  shift-invariant  if  /  €  S  implies  /(•  +  j)  is  in  S  for  every  multi-integer  j  £  ld.  These 
works  characterize  the  structure  of  shift-invariant  spaces  and  utilize  this  structure  to  settle 
several  questions  in  approximation  and  the  construction  of  wavelets. 


It  is  shown  for  example  how  the  closed  shift-invariant  subspaces  of  Z,2(Rd)  can  be  charac¬ 
terized  by  their  Fourier  transforms  and  how  to  describe  the  L 2(Rd)  projectors  onto  finitely 
generated  shift -invariant  spaces.  This  characterization  can  then  be  used  to  solve  several 
other  problems.  For  example,  for  an  arbitrary  closed  shift  invariant  subspace  S  of  Z^R**), 
it  is  characterized  when  the  dilates  Sh  of  S  provide  Z.2(Rd)-approximation  order  0(hk)  for 
all  functions  functions  /  in  the  Sobolev  space  \Vk(L2(Rd)).  This  extends  and  completes 
earlier  investigations  of  this  type  of  approximation  inititiated  by  Strang  and  Fix,  in  that 
no  assumptions  are  made  on  S.  Applications  to  wavelets  extend  and  complete  earlier  in¬ 
vestigations  of  Micchelli,  Chui  and  Wang,  Jia  and  Micchelli,  Riemenschneider  and  Shen, 
Lorentz  and  Madych,  Chui,  Stockier,  and  Ward,  and  others.  The  main  contribution  is  that 
all  a  priori  conditions  of  growth  and  smoothness  on  the  generating  functions  in  wavelet 
constructions  can  be  eliminated. 

There  are  several  obstacles  to  the  application  of  wavelets  for  problems  on  a  finite  do¬ 
main  caused  by  the  spill  over  (outside  the  domain)  of  the  supports  of  these  wavelets.  This 
is  one  of  the  obstacles  in  applying  wavelets  to  partial  differential  equations  and  to  the 
use  of  wavelets  in  segmentation  procedures  for  image  processing.  Recent  work  of  the  Pi’s 
and  their  collaborators  have  given  wavelet  construction  which  eliminate  some  of  these  dif¬ 
ficulties.  In  [28]  A.  Cohen,  I.  Daubechies,  B.  Jawerth,  and  P.  Vial  give  a  construction 
of  wavelets  on  the  unit  interval  [0.1].  The  construction  is  very  simple  and  robust,  and 
leads  to  fast  algorithms.  In  [2]  L  Andersson,  N.  Hall,  B.  Jawerth,  and  G.  Peters  extend 
these  results  in  different  ways.  W<*  consider  arbitrary  intervals,  which  is  important  for 
image  analysis  and  PDEs,  and  recursive  generation  of  wavelets  through  general  orthonor¬ 
mal  sets.  We  also  show  how  these  constructions  lead  to  smoothness  preserving  extension 
theorems,  important  for  PDE's  and  CAD,  and  how  they  can  be  used  to  construct  wavelets 
satisfying  boundary  conditions.  The  latter  has  been  an  open  problem  for  some  time.  In 
[5],  Andersson,  Dahlberg,  and  Jawerth  give  similar  constructions  for  general  closed  sets 
(including  sets  with  Lipschitz  boundary).  The  closed  set  wavelets  and  related  ideas  are 
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ideal  for  parallel  processing  and  this  is  exploited  by  Huntsberger,  Jawerth,  Metlapalli  and 
Peters  in  [22]. 

Our  study  of  nonlinear  multivariate  wave  propogation  [12]  described  below  has  pointed 
to  several  deficiencies  in  the  usual  construction  and  application  of  wavelets  in  several  space 
dimensions.  In  particular  we  observed  that  the  usual  wavelet  bases  are  deficient  in  being 
able  to  resolve  multivariate  shocks  because  their  dyadic  scaling  is  the  same  in  each  coor¬ 
dinate  direction.  To  circumvent  these  drawbacks,  we  have  introduced  hyperbolic  wavelets 
in  [26]  which  allow  for  variable  scaling  in  coordinate  directions.  The  theoretical  properties 
of  these  nonisotropic  wavelet  bases  were  developed  in  [26]  including  their  use  in  linear 
approximation  and  nonlinear  compression.  We  are  currently  investigating  application  of 
these  bases  to  image  processing. 

While  the  usual  wavelet  basis  and  the  hyperbolic  basis  are  two  contrasting  alternatives 
for  wavelet  decompositions,  there  is  actually  a  plethora  of  wavelet  bases  which  mitigate 
between  these  two  extremes.  In  this  direction,  we  have  investigated  in  [27],  the  adaptive 
selection  of  bases  with  variable  coordinate  scaling  for  applications  to  image  processing 
described  below. 

Partial  Differential  Equations.  We  have  concentrated  on  two  of  the  most  prominent 
problems  in  partial  differential  equations:  nonlinear  equations  such  as  conservation  laws 
which  arise  in  the  study  of  wave  propogation  and  elliptic  problems  on  nonsmooth  domains 
which  arise  in  scattering  theory  and  many  aspects  of  computer  aided  design. 

Ron  DeVore  and  Brad  Lucier  have  studied  regularity  theorems  and  numerical  methods 
for  solving  nonlinear  conservation  laws.  The  solution  to  such  equations  develop  discontinu¬ 
ities  (so  called  “shocks”)  even  for  C°°  initial  conditions.  Their  earlier  results  showed  that 
the  Besov  spaces  B“(Lr),  r  :=  (a  +  l)-1  are  regularity  spaces  for  scalar  conservation  laws 
in  one  space  dimension.  In  [12],  they  characterize  all  Besov  spaces  which  are  regularity 
spaces  for  such  conservation  laws.  They  consist  of  the  above  Besov  spaces  and  certain 
other  Besov  spaces  with  a  <  1.  These  results  can  also  be  appliedto  multivariate  scalar 


conservation  laws.  They  show  that  none  of  the  Besov  spaces,  with  a  >  1,  can  be  regularity 
spaces  for  conservation  laws  in  severed  space  dimensions.  The  importance  of  these  results 
is  that  it  indicates  that  conservation  laws  in  several  space  dimensions  must  involve  non- 
classical  smoothness  spaces.  An  indication  of  the  structure  of  these  non-classical  spaces  is 
also  given. 

Dahlberg  and  Jawerth  and  their  collaborators  have  also  studied  numerical  methods  for 
the  fast  solution  of  elliptic  PDE’s  on  domains.  In  [29]  we  have  established  that  the  Galerkin 
method  can  be  used  for  solving  integral  equations  corresponding  to  classical  elliptic  PDE’s 
on  nonsmooth  domains.  In  fact,  we  have  shown  that  the  mesh  size  can  be  chosen  in  an 
adaptive  way  controlled  by  local  curvature  properties  of  certain  approximate  domains.  We 
have  also  constructed  finite  element  spaces  based  on  adaptive  constructions  of  wavelets. 
These  constructions  were  precursors  to  those  that  lead  to  wavelets  on  closed  sets  and  can 
now  be  improved. 

Image  processing.  Ron  DeVore,  Bjorn  Jawerth,  and  Brad  Lucier  [13],  [14]  have  inves¬ 
tigated  algorithms  for  image  and  surface  compression  based  on  wavelet  decompositions. 
The  algorithms  developed  in  these  papers  have  been  shown  by  R.  DeVore,  G.  Kyriazis, 
D.  Leviatan,  and  V.  Tikhomirov  [16]  to  be  optimal  methods  for  compression  in  the  sense 
of  widths  or  optimal  recovery.  R.  DeVore  and  B.  Lucier  [17]  have  derived  new  algorithms 
for  noise  removal  based  on  wavelet  decompositions.  The  algorithms  for  noise  removal  can 
be  shown  to  be  optimal  in  the  sense  of  nondeterministic  optimal  recovery.  These  results 
compliment  parallel  investigations  of  Dnnoho  and  Johnstone  on  wavelet  shrinkage. 

In  [23],  L.  Hsiao,  B.  Jawerth.  B  Lucier,  and  X.  Yu  have  shown  optimal  compression 
algorithms  for  spaces  with  a  Littlewood  Paley  description,  such  as  Lp ,  Sobolev  spaces, 
BMO ,  certain  Besov  spaces,  etc.  These  results  improve  and  extend  some  in  [14].  In  this 
paper,  optimal  algorithms  for  compression  of  matrices  are  also  established.  These  results 
have  been  implemented  numerically  by  two  graduate  students  and  a  systematic  comparison 
with  the  results  of  Beylkin-Coifman-Rokhlin,  The  fast  wavelet  transform ,  has  been  carried 
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out  (the  results  are  now  part  of  the  students'  theses). 

In  [25]  B.  Jawerth,  A.  Laine,  and  A.  Tirumalai  have  developed  a  fast  technique  for 
magnifying  and  reducing  images  based  on  wavelet  methods.  In  contrast  to  the  standard 
method,  so  called  pixel  replication,  this  method  does  not  introduce  edge  artifacts.  In 
addition,  it  commutes  with  wavelet  compression  in  the  sense  that  the  magnification  can 
be  done  first  and  then  the  compression  or  vice  versa. 

We  are  developing  fast  algorithms  in  image  processing  for  their  application  to  digital 
mammography  and  the  early  detection  of  breast  cancer.  These  include  fast  algorithms  for 
compression  of  mammographic  images  [27]  which  do  not  lose  any  of  the  essential  features 
of  the  mammogram;  fast  methods  for  noise  removal;  and  adaptive  algorithms  for  feature 
classification  including  microcalcification  clusters,  inarginated  masses,  and  arcliitectural 
distortions. 

Computer  aided  design  We  have  continued  our  development  of  the  second  generation 
of  the  CAD  system  SLIP.  We  now  have  a  prototype  running  which  allows  us  to  reconstruct 
a  3D  CAD  description  from  photographic  images  of  a  manufactured  part. 

2.  Collaboration  with  industry.  We  are  continuing  our  collaboration  with 
various  industries,  government  laboratories,  and  university  centers. 

We  are  continuing  our  close  collaboration  with  Volvo,  Sweden.  This  collaboration  in¬ 
volves  CAD,  optimal  design  and  tool  design.  We  are  also  working  with  Metimur  and  Volvo 
concerning  computer  representation,  and  manufacturing,  of  industrial  objects  from  their 
photographic  images.  One  of  tin*  goals  with  this  project  is  to  allow  the  use  of  photo¬ 
graphic  images  for  quality  control  and  verification  of  specifications  of  parts  supplied  by 
subcontractors. 

A  collaboration  with  several  of  the  companies  within  McDonnell-Douglas,  St.  Louis, 
has  also  started.  The  collaboration  involves  PDE  solvers,  image  compression,  ATR,  and 
transient  analysis.  The  companies  involved  are  McAir,  MD  Missile  SC  ,  and  MD  Research 
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Laboratories. 

Collaboration  is  also  ongoing  with  Magnavox  on  video  compression  and  corresponding 
hardware  design. 

We  have  also  had  initial  discussions  with  United  Technologies,  concerning  PDE  solvers, 
manufacturing  and  CAD,  and  with  KSR  concerning  video  compression. 

We  collaborate  with  Hercules  Corporation  on  the  develoment  of  wavelet  based  algorithms 
for  Automated  Target  Detection.  This  research  has  centered  around  preprocessing  using 
wavelet  decompositions  to  be  followed  by  a  nonlinear  neural  network  for  the  detection  and 
classification. 

With  the  University  of  South  Florida  Cancer  Institute  we  collaborate  on  the  applications 
of  our  fast  image  processing  algorithms  in  medical  imagery  including  the  early  detection 
of  breast  cancer  through  digitized  mammography. 

Our  close  collaboration  with  SMTC  is  also  continuing.  SMTC  is  a  national  center 
at  USC.  Its  objective  is  to  provide  small  and  midsize  manufacturing  companies  in  the 
Southeastern  region  with  access  to  start  of  the  art  CAD/CAM  technology.  Currently, 
about  400  companies  are  associated  with  SMTC.  Our  work  with  SMTC  concerns  wavelet 
compression  methods,  fracture  mechanics,  CAD,  and  quality  control  through  contact  free 
measuring  techniques. 

The  activities  which  were  initiated  by  this  AFOSR  grant  are  also  continuing  as  part  of 
a  project  funded  by  a  new  DARPA  grant  for  Surface  modelling  and  Fast  PDE  solvers  (Pis 
B.  Dahlberg  and  B.  Jawerth).  This  new  project  involves  SMTC  as  well  as  the  Computer 
Science  and  Mechanical  Engineering  Departments. 

The  Industrial  Mathematics  Initiative.  Within  the  Department  of  Mathematics  at 
USC,  we  have  established  a  semi-autonomous  unit,  "The  Industrial  Mathematics  Initia¬ 
tive”  (IMI).  The  goal  of  IMI  is  to  enact  the  rapid  transfer  of  advanced  mathematics  and 
numerical  analysis  into  the  industrial  sector.  We  have  recently  also  been  awarded  an  NSF 
EPSCoR  grant  in  support  of  IMI  for  projects  related  to  the  work  on  this  AFOSR  contract. 
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